Linear Stability of Dilatonic Black Holes* 



Panagiota Kanti 

Division of Theoretical Physics, Physics Department, 
University of Ioannina, Ioannina GR 451 10, Greece. 



Abstract 

In this talk, we recall the most important features of the Dilatonic Black Holes 
which arise in the framework of the Einstein-Dilaton-Gauss-Bonnet theory and which 
are dressed with a classical long range dilaton field in contradiction with the existing 
"no-hair" theorems of the Theory of General Relativity. We demonstrate linear 
stability of these black hole solutions under small spacetime-dependent perturbations 
by making use of a semi-analytic method based on the Fubini-Sturm's theorem. As 
a result, the Dilatonic Black Holes constitute one of the very few examples of stable 
black hole solutions with non-trivial "hair" that arise in the framework of a more 
generalised theory of gravity. 
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1 Introduction 



The Dilatonic Black Holes arise in the framework of the one-loop corrected four- 
dimensional effective theory of the heterotic superstrings at low energies. In ref. [|TJ, 
we have demonstrated numerically the existence of the Dilatonic Black Hole solutions 
with a regular event horizon and an asymptotically flat behaviour at infinity in the 
presence of higher-derivative gravitational terms such as the Gauss Bonnet (GB) 
curvature-squared term. According to our analytic arguments |l] , it was the presence 
of this term which led to the by-passing of the "no-scalar-hair" theorems |2j and the 
existence of black hole solutions dressed with non-trivial classical dilaton hair. 

The existence itself of the Dilatonic Black Holes is beyond any doubt since 
our results were rederived and confirmed by other research groups in subsequent 
works |§ {§. But, an important question related to the nature and fate of our 
black hole solutions arises next : the question of stability. Although the "no-hair" 
theorems refer to the existence of black hole solutions and not to their stability, 
there is a tendency in the literature to reject some of these solutions due to their 
instability under spacetime perturbations. At the same time, the conception that 
the "no-hair" theorems are indeed valid once we have assured the stability of the 
black hole solutions has started to form among the scientists. For these reasons, we 
need also to check the stability of our solutions, that is the Dilatonic Black Hole 
solutions 0. 

In the next paragraph, we are going to recall the most important features of 
the Dilatonic Black Holes while in paragraph 3 we shall present our semi-analytic 
method for the study of our solutions under linear spacetime-dependent perturba- 
tions H . The key feature of our method is the reduction of the time-dependent 
equations of motion of the theory to a single one-dimensional Schrodinger type dif- 
ferential equation. As we shall see, the absence of bound states in the spectrum 
of this equation corresponds to the absence of unstable modes or equivalent to the 
stability of our classical black hole solutions of the time-independent equations of 
the theory |jj . This linear stability, if it comes out to be true, will be an outstand- 
ing result since the corresponding stable black hole solutions will be unique in the 
framework of a generalised theory of pure gravity. 

2 Dilatonic Black Holes 

The action functional of the Einstein-Dilaton-Gauss-Bonnet (EDGB) theory, in 
the context of which the Dilatonic Black Holes have arisen, has the following form H 
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where 

k% b = R^R^ - m^nr + r 2 (2) 

is the well-known curvature squared Gauss-Bonnet term and stands for the dilaton 
field. For simplicity, we have ignored all the other scalar fields of the effective 
superstring theory, the axions a and b and the modulus field a, as well as the gauge 
fields. 

The dilaton field and the Einstein's equations which follow from the action ([1]) 
can be written in a covariant form in the following way 

-^=d^/=g^) = -^e*I^ B (3) 
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We assume that the spacetime background is spherically symmetric and it is 
described by the following line element 

ds 2 = e r ^dt 2 - e A ^dr 2 - r 2 (d6 2 + sin 2 9 dip 2 ) . (7) 

If the above ansatz is to describe a black hole solution, we must demand the 
following behaviour for the components of the metric tensor and the scalar filed 
near the horizon of the black hole 

e- A(r) = Xi(r-r h )+X 2 {r-r h ) 2 + ... 

e r(r) = lx (r- Th ) +72 ( r _ rfc )2 + ... ( 8 ) 

(j){r) = h + <j)' h (r-r h ) + ( j)'; i (r-r h ) 2 + ... 



Substituting the above expansions in the equations of motion, we find that the 
derivative of the dilaton field at the horizon satisfies the relation 
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which results in the following constraint for the coupling function of the dilaton field 
with the Gauss-Bonnet term 



< -7= • (10) 

If the coupling function does not satisfy the above constraint, the solutions of the 
equations of motion of the theory can no longer be described by the concept of the 
black hole. At the other end of the radial space, that is in the limit r — > oo, one 
assumes the following asymptotically flat behaviour: 
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The constants M and D that appear in the above expressions stand for the ADM 
mass and the dilaton charge of the black hole, respectively. 

If we integrate numerically the equations of motion by making use of 

the ansatz (0) and the expansions (^) and flTT|), we are led to a continuous one- 
parameter family of regular black hole solutions. The solution for the scalar field 
and the metric components is displayed in Figures 1 and 2. The main feature of 
the solution, that we are going to make use of during the stability analysis, is the 
monotonic, non-intersecting behaviour of the dilaton field from to infinity. 
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Figure 1 : Dilaton configurations for a Figure 2 : Metric components gu and 
family of black hole solutions, for g rr for the 4>h = — 1 and = 1 

fixed r/j = 1 and various values of (ph- black hole solution. 
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The existence itself of these black hole solutions, of the Dilatonic Black Holes, are 
in disagreement with the well-known "no-hair" theorems of the Theory of General 
Relativity. These theorems state that the only parameters that may characterize a 
black hole are the mass M, the charge Q and the angular momentum J and that the 
only long range fields that can be associated with the black hole are the gravitational 
and the electromagnetic ones. The Dilatonic Black Holes are clearly characterized 
by another long range field, that of the dilaton field, a feature which is prohibited 
by the Theory of General Relativity. As we mentioned before, it is the presence of 
the Gauss-Bonnet term that causes the by-passing of the "no-hair" theorems and 
the existence of the long range scalar field. On the other hand, this disagreement 
with the "no-hair" theorems is only partial in the sense that we have not found any 
new conserved quantity apart from the aforementioned ones. The dilaton charge D 
has turned out to be a dependent on the mass M quantity, thereby leading to the 
secondary nature of the dilaton "hair" ||. 



3 Linear Stability Analysis 

In this section, we are going to present our semi-analytic method for the study 
of our solutions under linear spacetime-dependent perturbations ||. For the needs 
of our analysis, we are going to assume that both of the metric components as well 
as the scalar field depend not only on the radial coordinate r but also on the time 
coordinate t. Moreover, according to the linear stability analysis |7| [[5], we may 
write the three unknown functions of the problem as a sum of two parts in the 
following way 

r(r,t) = r(r) + 5T(r,t) = T(r) + 5T(r)e iat 

A(r,t) = A(r) + 8A(r,t) = A(r) + 5A{r)e iat (12) 

0(r, t) = (f)(r) + 8(f>(r, t) = 0(r) + 5(f)(r)e iat 

The first part of each of the above functions is a purely radial part which corre- 
sponds to the classical Dilatonic Black Hole solutions already found in Ref . [[[]] . The 
second part is the product of a small radial part, denoted by the prefactor 5, by an 
harmonic function of time. These small radial parts are the linear perturbations of 
our solutions. Note the appearance of the parameter a in the above expressions. It 
is easy to understand that for negative a 2 , that is for purely imaginary a, even small 
perturbations may be led to extremely large values at late times. For this reason, 
we have to assure the absence of negative a 2 in the case of our black hole solutions. 

If we substitute the expressions in the time-dependent equations of motion 
(ED - (ED an d make use of the harmonic dependence of the perturbations on time t, we 
end up with a group of differential equations with respect to the radial coordinate 
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r for the linear perturbations 5T, 5 A and 5(f). After some long and tedious algebraic 
computation, we manage to construct a single differential equation for the dilaton 
perturbation 5(f> from which the other two perturbations, 5T and 5 A, have been 
eliminated. This differential equation has the following structure 

A5(f)" + 2B 5(f)' + C 5(f) + a 2 E 5(f) = (13) 

where A, B, C, and E are rather complicated functions of (f>, 0', 0", A, A', T, V and 
r". All of these coefficients are very well behaved near infinity, that is in the limit 
r — > oo, according to the following expressions 

A = 1 + 0(± E ) , B = l + ^ + 0(±) (14) 
D 2 1 AM AM 2 / 1 \ 

c =^°^) > . as) 

However, three of them, B, C and i£, take on infinite values near the horizon of the 
black hole, while A remains finite, in the following way 

* - TT§ + °(— ) ■ ^(^ + 0(1) (i6) 

- V , + 0(+.) (17) 
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where x = 1 — 6e 2r ^ h /r^. 

As a result, the differential equation ( |13"D is not well defined at one point only of 
the radial space, that is at the event horizon r^. In order to remove this singularity, it 
is necessary to introduce a new radial coordinate, the so called "tortoise" coordinate 
r*, which is related to the old one in the following way 

±1 = e -(r-A)/2 (19) 

The introduction of the new radial coordinate leads to the extension of the radial 
space [r/i, oo) over the entire real axis (— oo, oo) in such a way that the relative slope 
of any two curves in Fig. 1 remains unchanged. Nevertheless, the introduction of the 
new coordinate is powerful enough to render all of the coefficients in the differential 
equation fll3"D finite over the entire radial space. Now, the perturbed equation for 
the dilaton field takes the form 

^ + 2B^ + (C + ^)^0 (20) 
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where 

A = A , B = Be^ 2 -^ d{T ~ h) , C = e r ~ A C , S = e r ~ A E (21) 

4 dr* v ' 

For reasons that will become clear later, we have to make another step in order to 
eliminate the term in equation (|20D which is proportional to 5(f)' . This step involves 
the use of the auxiliary function F defined as 




exp / -rdr*'\ (22) 



and the definition of a new function u as the product of the auxiliary function by 
the dilaton perturbation, u = F 5<p. Then, the differential equation for 8(f) takes the 
form 



C 2 £ B 2 (B 
A + a A-A*- p * {A 



u = (23) 



where 



P* = A • C 24 ) 
dr* 



In this form, equation (P3|) is an ordinary Schrodinger equation with well defined 
coefficients over the whole radial space. At the same time, this differential equation 
is an eigenvalue problem where different values of the parameter a 2 correspond to 
different eigenvalues. It is easy to understand that, since the number of negative 
values of a 2 is equal to the number of the unstable modes of our black hole solutions, 
the absence of states with negative a 2 in the above eigenvalue problem is equivalent 
to the linear stability of the Dilatonic Black Holes. 

In order to demonstrate the absence of the aforementioned states from the spec- 
trum of the eigenvalue problem (|2"3"|), we are going to make use of the asymptotic 
behaviour of the eigenfunctions u at both limits of the radial space, r — > and 
r — > oo, and of the Fubini-Sturm's theorem ||. We are going to concentrate our at- 
tention on two different kinds of eigenfunctions, namely on m which corresponds to 
zero eigenvalue a 2 = and on u a which corresponds to negative eigenvalues a 2 < 0. 
We start by noting that near the horizon of the black hole, that is in the limit r — > Vh 
or equivalently r* — > — oo, the differential equation ( ^3|) assumes the form 



V \u a + k 2 u a = (25) 



k 2 = 1 +a 2 = k 2 + a 2 (26) 

(1 + Vl-6e 2 **)Vl-6e 2 *" v ; 



where 



and where we have made use of the definitions (pip and the asymptotic expansions 
(§). It can be easily seen that, while for a 2 = the corresponding eigenfunction uq 
takes on a constant, non-zero value at the event horizon, for —k% < o 2 < the only 
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acceptable value of u a , in order to ensure the finiteness of 5(p, is : «„ ~ e' fe ' r * — > 0. 
Here, we have to add that the eigenfunctions with a 2 < —k 2 < are absent from 
the spectrum of the states since they have been shown to violate the continuous, 
non-degenerate nature of the unbound states ||. 

At the other end of the radial space, in the limit r — > oo, it can be easily seen 
from the expressions ([T^)-([T5|) that the differential equation (|23|) takes the simple 
form 

plu a + a 2 u a = . (27) 

In the same way, the only acceptable behaviour of the eigenfunction u for a 2 < 0, in 
the limit r — > oo, is : u a ~ e~' cr ' r ' — > while for a 2 = 0, the solution is Uq = C\ r + c 2 - 
Note that only the eigenfunctions u a that correspond to negative values of the 
parameter a 2 vanish at both limits of the radial space. For this reason, only these 
eigenfunctions correspond to physical perturbations of our black hole solutions |7| . 



Now, we arrive at the final step of our stability analysis. This final step involves the 
use of the Fubini-Sturm's theorem of ordinary differential equations ||. According 
to the theorem, we consider the following two differential equations : 

< + (gi-p?-pi)«i = (28) 

«2 + (?2-^-j/ 2 )«2 = (29) 
If the coefficients of the above equations satisfy the following relation 

p 2 + p\ ~ Q2 < p[ +p\-qi (30) 

throughout the interval [a, b], then, between any two consecutive zeroes of function 
Ui, in the interval [a, b], there is at least one zero of function U2- 

If we apply the above theorem to the spectrum of bound states of a Schrodinger 
type eigenvalue problem, we obtain the well-known "node rule" , according to which, 
if we arrange all the bound states in an increasing order of eigenvalues, the n - th 
eigenfunction has n — 1 nodes ||10|| . This means that the ground state in the spectrum 



of the bound states, denoted by u^, has no zeroes throughout the interval [a, b]. 



Finally, we consider the case where u\ = Ub with a 2 being the most negative 
eigenvalue of the system and u 2 = uq with a 2 = 0. We multiply equation ( |2"8"|) by u 2 
and equation (|29|) by u\, subtract the two equations and integrate the result over 
the entire domain of r*. Then, we obtain : 



[u a p*u - u p*u b/ 



- a 2 ) I ^-u b u Q dr* (31) 



By making use of the asymptotic behaviour of the eigenfunctions uq and u b at 
both limits of the radial space, a behaviour which was discussed previously, it can 
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Figure 3 : The coefficient S/A for the = 1 and 0^ = — 1 black hole solution. 
The definite sign of this coefficient is obvious. Also shown is the coefficient E/A 
before the introduction of the "tortoise" coordinate. 



be easily seen that the left hand side of the above equation is zero. The same 
must also hold for the right hand side of this equation. Note that the integral that 
appears on the right hand side contains the product of three functions, S/A, Ub and 
Uq. The ratio of the coefficients S and A, after the introduction of the "tortoise" 
coordinate, is a well defined function over the entire radial space as it is shown in 
Figure 3. Moreover, it is a positive definite function without any zeroes throughout 
the radial domain. On the other hand, according to the "node rule" , the function Ub 
is the ground state in the spectrum of bound states of the eigenvalue problem and, 
as such, it has no zeroes as well. The function w demands some more attention 
since it corresponds to a vanishing value of the parameter a. In this case, the time 
dependence in the expressions (0) of the linear perturbations disappears and the 
small radial parts can be absorbed into the large radial parts which describe the 
classical Dilatonic Black Hole solutions. As a result, the eigenfunction uq can be 
easily constructed out of the difference of any two of the curves in Figure 1. The 
monotonic, non-intersecting behaviour of these curves ensures the absence of zeroes 
in the expression of Uq. Since each one of the three functions that appear inside 
the integral has a definite sign, their product will have a definite sign as well and 
the integral can never vanish. Then, the only consistent case is the degenerate one 
a 2 = Uq = 0, which means that, actually, the ground state of the spectrum is the 
one with the zero eigenvalue and that all of the bound states with negative values 
of the parameter a 2 are absent. Since the number of negative values of a 2 is equal 
to the number of the unstable modes of our solutions, the absence of bound states 
implies the linear stability of the Dilatonic Black Holes. 
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4 Conclusions 



By the use of a semi-analytic method, we have demonstrated the linear stability 
of the Dilatonic Black Holes, which arise in the framework of the Einstein-Dilaton- 
Gauss-Bonnet theory, under small spacetime-dependent perturbations. As a result, 
our black hole solutions have no reason to be rejected, even if we make use of 
the argument of stability, and the concept of the validity of the "no-hair" theorem 
among the stable solutions has proven to be false. Our result is extremely important 
since it renders our black hole solutions linearly stable and makes the long range 
dilaton field, or equivalently the "dilaton hair", one of the very few examples of 
stable, although of "secondary type" , hair that bypasses the "no-hair" theorems of 
the Theory of General Relativity. 

The stability itself of the Dilatonic Black Holes under spacetime-dependent per- 
turbations is very important since the lifetime of these black hole solutions becomes 
substantially longer. As a result, Dilatonic Black Holes, that may have been formed 
in the past, are rather possible to have survived until the present epoch. In this 
case, the possibility of detecting them is not negligible at all and this might provide 
a test for the predictions of superstring theory at low energies. 

As a final remark, we would like to note that the Schwarzschild Black Hole, 
which arises in the framework of the Theory of General Relativity, is also linearly 
stable. This means that our black hole solutions and the Schwarzschild Black Hole 
share the same kind of stability. Moreover, the Dilatonic Black Holes arise in the 
framework of the Einstein-Dilaton-Gauss-Bonnet theory which is a pure gravity 
theory-the dilaton field simply plays the role of the connecting link between the 
scalar curvature R and the Gauss-Bonnet term. Taking into account the above as 
well as the fact that these two black hole solutions are characterized by the same 
parameter at infinity, the mass M, we are led to the conclusion that the Dilatonic 
Black Holes can be considered as the most direct generalisation of the Schwarzschild 
Black Hole in the framework of the effective superstring theory. 
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